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ABSTRACT 


W 


Forced convection and subcooled boiling heat transfer data for liquid 
nitrogen and liquid neon were obtained in support of a design study for a 
30 tesla cryomagnet cooled by forced convection of liquid neon. This de- 
sign precludes nucleate boiling in the flow channels as they are too small 
to handle vapor flow. Consequently, it was necessary to determine boil- 
ing incipience under the operating conditions of the magnet system. 

The cryogen data obtained over a range of system pressures, fluid 
flow rates, and applied heat fluxes were used to develop correlations for 
predicting boiling incipience and convective boiling heat transfer coeffi- 
cients in uniformly heated flow channels. The accuracy of the correlating 
equations was then evaluated. 

A technique was also developed to calculate the position of boiling 
incipience in a uniformly heated flow channel. Comparisons made with 
the experimental data showed a prediction accuracy of ±15 percent. 


INTRODUCTION 

The data reported herein were obtained in support of a design study 
for a 30 tesla cryomagnet described in reference 1. Because the coolant 
channels of the magnet coils are too small to handle vapor flow, it was 
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necessary to determine boiling incipience and convective boiling coeffi- 
cients under the operating conditions of the magnet system. Forced con- 
vection heat transfer data obtained for this design were correlated and are 
presented in reference 2. 

In spite of the voluminous literature on boiling heat transfer, the pre- 
dictability of boiling incipience and convective boiling in flowing heat trans- 
fer systems is still quite limited. Available empirical and semi-empirical 
techniques (presented in refs, 3, 4, 5) are mostly derived from water data. 
In this report we extend these data to liquid neon and liquid nitrogen. 

The convective boiling data were correlated by a technique similar to 
that of reference 5. In addition, the position of boiling incipience in the 
heated channel was calculated from the value of the correlation parameter 
at the incipient boiling point. Comparisons were then made with the exper- 
imental data. 


NOMENCLATURE 

A area, cm^ 

Cp specific heat, J/g-K 

Cito.Ce constants in eq. (1) 
d diameter , cm 

G Pb v b = W / A 

H enthalpy, J/g 

2 

h heat transfer coefficient, W/cm -K 

k thermal conductivity, J/sec-cm-K 

L distance along test section, cm 

Nu Nusselt number, hd/k 

Pr Prandtl number, C jtx/k 

pressure, MPa bars 


P 
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Q total heat input, W 
q heat flux, W/ciu 
Re Reynolds number, pVd/ jj, 

T temperature, K 
V velocity, cm/sec 

x^ heated length for boiling incipience, cm 
X heat of vaporization, J/g 
l± viscosity, g/cm-sec 

o 

p density, g/cm 
w mass flow rate, g/sec 
Subscripts 
b bulk fluid 
boil boiling- 
calc calculated 
fc forced convection 
ht heat transfer 
in inlet 
1 liquid 
s saturated 
v vapor 
w wall 


APPARATUS AND PROCEDURE 

The entire heat transfer apparatus, composed of the cryogen tank, test 
section, and radiation shield, was housed in a steel vacuum tank (fig. 1), 
The instrumented test section, designed to simulate the small flow channels 
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of the proposed liquid neon cooled magnet, was installed to permit vertically 
upward flow and was uniformly joule heated. Pressurization of a cryogenic 
supply tank or "blowdown” was used to force the fluid to flow through the 
test section. 

The joule heated Inconel test section was 25 cm long by 0. 198 cm i.d. 
and 0. 256 o, d. with a resistance of about 0. 148 £2. An ac power supply pro- 
vided up to 6. 8 V across the test section for a maximum heat flux of 21.2 
2 

W/cm . Tube wall temperature measurements were made with 40 gage 
Chromel-constantan thermocouples located at positions 5. 7, 8. 2, 10. 7, 

13. 3, 15. 7, 18.2, 20. 7, and 23. 1 cm from the test-section inlet-side elec- 
trode. The thermocouples were referenced to the fluid temperature in the 
inlet mixing chamber. Fluid temperatures at the entrance and exit of the 
test section and in the cryogen tank were measiued using platinum resistance 
thermometers. Mixing chamber pressure taps were provided for pressure 
measurements. The flow rate was determined by using an orifice in con- 
junction with measurements of local temperature and pressure near the bot- 
tom of the dip tube in the cryogen tank as shown schematically in figure 1. 

Before taking any cryogenic data, the flow system and instrumentation 
accuracy of the test rig were initially checked by taking gaseous nitrogen, 
neon, and helium heat transfer data. The forced convection heat transfer 
behavior of these gases is well known, so the experimental heat transfer co- 
efficients from these tests provided adequate information to judge the reli- 
ability of the instrumentation. Since these results were within ±10 percent of 
predictions by the Dittus-Boelter equation the system was considered 
operational. 

With liquid cryogens as the working fluid, the cooldown procedure was 
tailored to the type of fluid used. When the working fluid was liquid nitrogen, 
the flow system was initially cooled by the fluid itself. Liquid neon, on the 
other hand, is too costly to be used as a cooldown fluid. To minimize boil- 
off Kvises, the rig was precooled to approximately 20 K with cold helium gas 
IhyLOre any transfer of liquid neon was made.. 

The test procedure required control of system pressure, flow rate, and 
heat input. After the initial cooldown, the supply tank was filled with the 
cryogen and tlidU pressurized. The pressurizing gas, either helium or neon, 


! 
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was also precooled to approximately 80 K with liquid nitrogen to minimize 
heat transfer to the working fluid. The system pressure and flow rate were 
set by control of the inlet and outlet valves. At each power setting, time 
was allowed for the system to reach steady -state conditions before the data 
were recorded. The same procedure was repeated over a range of flow 
rates and system pressures. 

FORCED CONVECTIVE BOILING CORRELATION 
Papell (ref. 5) found the relationship 



to be effective in the correlation of subcooled boiling water data when the 
equation constants were C^ = 90, Cg = 1.0, Cg = 1.2, C^ = 1.08, and 
C 5 = 0. 7. Preliminary data have shown that for liquid nitrogen and liquid 
neon modification of this relationship is required. 


Liquid Nitrogen Data 

The liquid nitrogen data are presented in figure -2 in terms of the ratio 
of boiling to nonboiling heat transfer coefficients (h boil /hf C ) plotted as a 
function of the three dimensionless parameters of equation (1): 
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An analysis of the data showed a best fit was obtained when was changed 
to 1,4 with C 2 and C 3 remaining the same as for the water data. In addi- 
tion, the forced convection heat transfer coefficient (h fc ) in equation (1) was 
calculated from a Dittus-Boelter type correlation that has been verified in 
reference 3 for both liquid nitrogen and liquid neon as 

Nu fc = 0.023(Re)°' 8 (Pr) 0 - 4 (3) 

with fluid properties at bulk temperature conditions. 

The data in figure 2 cover a range of test conditions that include system 
pressures from 1. 11 to 2. 87 MPa (11. 1 to 28. 7 bars), flow rates from 2. 7 to 
6. 5 g/sec, heat flux from 4. 1 to 21. 2 W/cm , and subcooling temperature 
differentials up to 26. 5 K. 

Figure 2 shows that the data fall along a line with a slope of 0. 7, which 
is the same slope as the reference 5 water data correlation. An analysis of 
the scatter shows that 88 percent of the data lie within ±10 percent of the cor- 
relating line. The constant equal to 100 fits the data and the final corre- 
lation becomes 



The ratio is a direct indication of the enhanced heat transfer 

properties of boiling to nonboiling. On figure 2 this ratio remains at a value 

of unity for all the nonboiling data and is greater than unity when boiling exists. 

Equation (4) is limited to sub cooled boiling because the parameter (\/H g - H^) 

becomes undefined as the bulk fluid becomes saturated. Data with 1— degree K 

2 

of saturated bulk conditions were not used because of excessive scatter. 

A frequent problem associated with data correlation is the suppression of 
significant information. Correlations involving exponential terms tend to sup- 
press the scatter of the data employed in developing the correlation equations. 
Consequently, appreciable uncertainties emerge when the correlation is used 
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to predict the effect of varying a parameter in the correlation. To test the 
accuracy of the data correlation to predict h given q, d>, and bulk condi- 
tions, equation (4) was used to compare calculated and experimental values 
of heat transfer coefficient ratios. The results (presented in fig. 3) show 
that 82 percent of the data fall within a ±10 percent scatter. The calculation 
of heat transfer coefficients from q, cl>, and bulk conditions is therefore con- 
sidered reasonably valid, with minor suppression of information. 

On the other hand, the calculation of heat flux (q) suffers because of its 
high dependence on temperature difference. * In order to determine the sen- 
sitivity of q calc , equation (4) was reformulated as a function of AT, oj, and 
bulk conditions: 



The heat flux was then calculated using equation (5) and compared with the 
experimental data. The results, presented in figure 4, show a data scatter 
of ±30 percent. Thus, the sensitivity of the heat flux calculation in equation 
(4) is about three times greater than that of the heat transfer coefficient 
ratio calculation. 


Liquid Neon Data 

Figures 5, 6, and 7 treat the liquid neon data in the same manner as the 
liquid nitrogen data. Figure 5 shows the correlation of neon data over a range 
of pressures from 1. 12 to 1. 66 MPa (11. 2 to 16. 6 bars), flow rates from 4„ 0 
to 10. 8 g/sec, heat flux from 1. 8 to 20. 5 W/cm^, and sub cooling differentials 
up to 11. 6 K. An analysis shows that 82 percent of the data lies within 


By analogy, the sensitivity is quite similar to that in pool boiling where 
q cc AT 3 . 
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±15 percent of correlation equation (4). The suppression of information (dis- 
cussed in the Liquid Nitrogen Data section) by this correlation, as determined 
by comparisons between calculated and experimental heat transfer coefficient 
ratios, is shown in figure 6, For a given set of q, co, and bulk conditions, 
approximately 75 percent of the neon data lie within a scatter of ±15 percent. 
The increased scatter in the neon data over the nitrogen data is partially due 
to the difference in temperature measurement accuracy. Temperature differ- 
ences used in the correlation for liquid neon were about half the values obtained, 
for nitrogen. In addition, less stable operation of the rig was noted when using 
neon rather than nitrogen. It is not understood, at this time, whether the in- 
stability was thermally or hydrodynamic ally driven. 

For a given AT, to, and bulk conditions, the calculated heat flux using 
equation (5) is compared •with experimental data in figure 7. An analysis shows 
that the data spread is about ±40 percent which again represents considerable 
suppression of information. 


BOILING INCIPIENCE 

The methods available in the literature to predict boiling incipience are 
described in references 3, 4, and 6. They require information such as the 
amount of superheat needed to initiate boiling or a knowledge of the bubble site 
cavity dimensions and the laminar sublayer thickness. All these parameters 
are functions of test conditions and are not readily available for most fluids. * 
The boiling correlation technique presented herein is premised on com- 
mercial tubes and can be used to determine local bulk enthalpy and. density 
conditions at boiling incipience for both liquid nitrogen and liquid neon. These 
local conditions can then be used to calculate the position of boiling incipience 
in the heated channel. 


* 

We do not mean to imply that sites, sublayers, etc. , are not important 
to the mechanism of nucleate boiling - only that for a given geometry they 
vary little with the process. 
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Correlation of Incipient Boiling Data 


An incipient boiling data correlation can be obtained directly from fig- 
ures 2 and 5 by inspection and definition. The noriboiling condition is repre- 
sented by a value of heat transfer coefficient ratio equal to unity and the boil- 
ing condition by the 0. 7 slope line. The intersection of these lines is defined 
as the incipient boiling point; for both cryogens it occurs at a value of the 
dimensionless parameters of 0. 0014. 

The incipient boiling correlation for both liquid nitrogen and liquid neon 
can be represented by the equation* 



For a given q, co, and system pressure level, to solve equation (6) for bulk 
enthalpy at the incipient boiling point requires an iteration technique since 
both bulk enthalpy Ha \ and bulk density p n \ are unknowns. 


L (bx.) 


^(bx^) 


Position of Boiling Incipience 

In a flow channel subjected to uniform convective heat fluxes, different- 
modes of heat transfer can exist simultaneously at different locations along 
the channel, and the transition position can be stable. For the transition from 
forced convection to nucleate boiling' the position of boiling incipience (x^) in a 
tube can be determined by the relation 


The apparent insensitivity of eq. (6) to temperature difference follows 
from the data. A given heat flux and pressure level fixes the wall temper- 
ature at boiling incipience. 
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x i = 


0 ) 


H. ^ - H. 
qhd V (bx i ) 


in ) 


(7) 


where equation (6) is used to calculate H^ x y 

The liquid nitrogen and liquid neon data Obtained for this study were used 
to calculate the position of boiling incipience over a wide range of test condi- 
tions. Experimental values of the position of boiling incipience were obtained 
from temperature profiles plotted as in figure 8; this figure shows the liquid 
nitrogen data at a system pressure of 2. 17 MPa (21. 7 bars) and a flow rate of 
6. 8 g/sec for three different heat fluxes. 

Whether nonboiling or boiling exists is evident from the slope of the lines 
drawn through the plotted points. The linear relationship between increasing 
wall temperature and heated length is typical of forced convection heat trans- 
fer data. Boiling occurs when the wall temperature is some finite value 
greater than the saturation temperature of the fluid which in figure 8 equals a 
superheat of 3 K. After boiling starts the wall temperature remains constant 
along the tube as long as nucleate boiling persists. 

The position of boiling incipience is located at the intersection of the 
forced convection and the nucleate boiling curves. For a heat flux of 21. 2 

O 

W/cm , figure 8 shows that boiling incipience occurs at a heated length of 
11 cm. The two lower heat flux curves are obviously nonboiling because of 
the monotonically increasing temperature profiles. The linear relationship 
between wall temperature and distance permits extrapolation to an intersec- 
tion with an extension of the boiling curve. In this manner the incipient boil- 
ing point for a uniform heat flux of 16. 8 W/cm with the geometry used in 
this study was found to equal 24. 6 cm of heated length and 48. 8 cm for a heat 
flux of 11, 4 W/cm^. 

The liquid nitrogen and liquid neon forced convection data and partial 
boiling data were treated in this manner. Comparisons were then made be- 
tween the calculated and measured incipient boiling points to determine the 
accuracy of the calculation procedure. The results are shown in figure 9 for 
both experimental and extrapolated data. All data fall within a ±15 percent 
scatter band. 
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SUMMARY OF RESULTS 

Liquid nitrogen and liquid neon forced convection subcooled boiling heat 
transfer data obtained for this study were used to develop correlations and 
techniques for predicting boiling incipience and convective boiling in uniformly 
heated flow channels. The results were obtained for the test conditions and 
flow geometry used in this study. 

The correlation equation for the subcooled boiling data showed that, given 
q, o>, and bulk conditions, 82 percent of the liquid nitrogen and liquid neon 
data correlated to within a scatter of ±10 and ±15 percent, respectively. This 
is considered a reasonably minor suppression of information in the calcula- 
tion of heat transfer coefficients. A comparison of calculated and experi- 
mental heat flux (q) for given AT, o>, and bulk conditions showed a scatter 
of ±30 percent for the nitrogen data and ±40 percent for the neon data. This 
represents considerable suppression of information, similar to that in pool 
boiling data. 

The correlation equation for boiling incipience, solved by iterating, was 
used to calculate the position of boiling incipience in a flow channel subject 
to uniform heat fluxes. Calculated values compared to experimental data and 
extrapolated values showed an accuracy to within ±15 percent. 
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Figure 1. • Liquid neon heat transfer facility. 
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Figure 2. - Liquid nitrogen convective boiling data corre- 
lation. Pressure ■ 1.11 to 2.8/ MPa (11.1 to 28.7 bars); 
flow rate ■ 2. 7 to 6. 5 g/sec; heat flux ■ 4. 1 to 21. 2 W/cm^. 
subcooling • 1.5 to 26.5 K. 



Figure 3. - Accuracy of liquid nitrogen data correlation ieq. (4i>. 
82 Percent of data within ±10 percent scatter. 




Figure 6. - Accuracy of liquid neon data correlation ieq. i4». 
75 Percent of data within +15 percent scatter. 



Figure 7. - Accuracy of correlation for heat flux up calcu- 
lation using equation i5i. Liquid neon data - 75 percent 
of data within ±40 percent scatter. 
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Figure 8. - Test section tube wall temperature profile data 
for liquid nitrogen at a system pressure of 2.17 MPa 
121.7 bars I and a flow rate of 6.8 g/sec. 
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Figure 9. - Accuracy of calculation procedure for predicting 
the position of boiling incipience Uj> in a uniformly heated 
tube using equations (61 and (7>. 
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